Quantum impurity models with the Coupled Cluster Method 
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We investigate tlie ground-state properties of the Anderson single impurity model (finite Coulomb 
impurity repulsion) with the Coupled Cluster Method. We consider different CCM reference states 
and approximation schemes and make comparison with exact Green's function results for the non- 
interacting model and with Brillouin-Wigner perturbation theory for the full interacting model. Our 
results show that coupled cluster techniques are well suited to quantum impurity problems. 
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The coupled cluster method (CCM) is a pow- 
erful method for investigating quantum many-body 
systemai— . Being derived from first principles, the CCM 
is universally applicable to many different fields, and 
its reputation, for being numerically accurate at rea- 
sonable computation costs, has been well established in 
nuclear physics^, quantum chemistry^i^, and quantum 
magnetism^i^i^. Although the method finds most appli- 
cation for discrete models, such as those of fields men- 
tioned above, it has also been succesfully applied to con- 
tinuum systems, such as in RefsJ^iii. In this paper, we 
consider the application of the CCM to the Anderson im- 
purity modcl^^, which consists of a single localised orbital 
coupled to a contiuum of electronic states. 



Since its inception, interest in Anderson-like impurity 
models has been high, as these models play important 
roles in strongly correlated systemsHiiiii^^ and trans- 
port through quantum dotsi^. Over the years, the An- 
derson model has been studied with many different ap- 
proachs: Bethe ansatz techniques are used to exactly 
describe ground-state propertie o^^i^^'^^ , and the Green's 
function method is used for an exact solution of the non- 
interacting casei^i^Sii^. Approximation methods, such as 
the variational method^^, projection operator o^^i^^ , Hub- 
bard operators^^^, and the numerical renormalization 
group method^! are also employed to the impurity mod- 
els; an approach similar to the CCM can as well be found 
in Ref,^. 



We will give a short introduction to the method, fol- 
lowed by an overview of the Green's function exact so- 
lution to the non-interacting Anderson model and self- 
consistent perturbation solution to the interacting one. 
We will then present a general disscusion of possible 
choices of CCM reference states and correlation opera- 
tors for the model. Several truncation schemes to these 
operators will be considered for both reference states and 
both models with no interaction and a finite Coulomb re- 
pulsion. The results will be compared with those of exact 
Green's function and self-consistent perturbation meth- 
ods. All of our calculations are done in equilibrium. 



I. COUPLED CLUSTER METHOD 

We begin with a brief outline of the CCM formalism; 
further details can be found in Refsi— . Consider a gen- 
eral many-body system described by a Hamiltonian H 
and exact ground-state eigenvector, j^*), such that 



H\^)=Eg\-^). 



(1) 



Our system may be described in terms of a reference state 
(or cyclic vector) |$) and a corresponding complete set of 
mutually commuting multiconfigurational creation oper- 
ators I C'j I . The set | C"] | is defined with respect to the 

reference state, such that C/ 1$) = = ($|C], V / 7^ 0, 
in a notation in which Cq = 1, the identity operator. 
In general, / is a set index and the operators C\ con- 
tain products of single-particle operators. The set in- 
dex {/} is complete in the sense that the set of states 

|c| |$)| provides a complete basis for the Hilbert (or 

Fock) space. The reference state, |(f>), must be chosen to 
be non-orthogonal to the actual wavefunction of the sys- 
tem, ($ l^*) 7^ and thus |<f>) plays the role of a vacuum 
state with respect to a suitable set of (mutually commut- 
ing) many-body creation operators {Cj}. For example, 
in describing an interacting Fermi gas, we might take the 
reference state to be the filled (non-interacting) Fermi 
sea, with the set of operators {Cj} creating all possible 
electron-hole excitations in this sea^^^. 

Within the single-reference CCM, the ket state of 
Eq. ([1]) is parametrised as 



|*)=e^|$) ; 5 = 5^.,Cl, 



(2) 



with S the CCM correlation operator and {s/} complex 
coefficients. This exponentiated form of the ground-state 
CCM parametrisation of Eq. ([2]) ensures the correct 
counting of the independently-excited correlated many- 
body clusters with respect to |<f>) that are present in the 
exact ground state \^). It also ensures the exact incorpo- 
ration of the Goldstone linked-cluster theorem^S^ which 
itself guarantees the size-extensivity of all relevant exten- 
sive physical quantities. 
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The eigen bra-state of our Hamiltonian is defined via 
the Schrodinger equation 



(vE-lff = Eg{^\. (3) 



Within the normal CCM^, this bra state is parame- 
terised independently as 



($|5e- 



5 = 1 



1^0 



(4) 



The exact groundstate eigcn-bra is given by the Hermi- 
tian adjoint of state |4>), but if the correlation opera- 
tor S is truncated, the approximate eigenket |$) may 
be no longer normalisable with itself and the adjoint- 
state ill-defined. We note that although Hermiticity 
for a truncated S is lost, the normalisation conditions 
(^flxj/) = (^Ivt/) = (<I>|<i>) = 1 are explicitly imposed. 

The ground-state properties of the system, then, are 
completely described by the set of CCM correlation coef- 
ficients {s/. Si} which are regarded as independent vari- 
ables. An arbitrary operator A will have a ground-state 
expectation value given as 

A EE = ($|5e-^Ae^|$) = A ({sj. Si}) . (5) 

In particular, the ground-state energy expectation func- 
tional H{{si, Si}) is given by 

ii ({s7. Si}) EE = ($|S'e"^i7c^|$). (6) 

By requiring _ff ({.s/. Si}) to be stationary with respect to 
variations in each of the (independent) correlation coef- 
ficients, one arrives at an expression for the ground-state 
energy 



Eg=Eg{{si}) = {<S>\c~''Hc^\<S>), 



(7) 



and the following coupled set of equations for the coeffi- 
cients 

SH/SSi=0 => (^IC/c-^i/C^I^) = 0, 1^0 ; (8) 
6H/Ssi=0 ($|S'e-^[iJ,C|]c^|$) =0, / 7^ .(9) 

The similarity transforms may be evaluated with help of 
the identity 



e-^He^ ^ H +[H,S] + [[H, S] , S] 



(10) 



which is guaranteed to terminate since H will contain 
only a finite number of annihilation operators C/ . 

It is important to notice that this (bi-)variational for- 
mulation does not provide an upper bound for Eg when 
the summations for S and S are truncated, due to the 
lack of exact Hermiticity when such approximations are 
made. 



II. ANDERSON MODEL 

We will study here the Anderson model of a single 
impurity coupled to a single reservoir. This model was 
originally introduced to describe magnetic impurities in 
metals^, and also finds application in describing trans- 
port through quantum dots^^. The Anderson Hamilto- 
nian reads 

k.cr a 

+ ^(Ffe4cfc. + F;cLd.), (11) 

where d\ is the creation operator of an electron of spin 
CT on the localised level and c^^ is the creation operator 
of a continuum electron with quantum numbers fc and cr. 
The energy of the dot level is e^, is the energy of a 
continuum level, TJ is the Coulomb interaction between 
two impurity electrons on the local level and T4 describes 
the coupling between dot and continuum state k. For 
simplicity, we have assumed that ea,-, Sk and Vk are spin- 
independent. 

We further specify the continuum as being a band of 
width 2W with centre chosen as our energy zero. We 
chose the number of electrons in the system, N^, to be 
that given by filling the band from to at zero 

temperaturci^ 

The important quantity that describes the hybridisa- 
tion of the localised level with the reservoir is the half 
tunnelling rate, 



r(e)^^.9(e)1^4)y,(,). 



(12) 



where g{e) is the density of states of the band, g{e) = 
J2k " ^k)- In obtaining numerical results, we will as- 
sume that this rate is constant across the energy range of 
interest, r(e) = F, although our technique is not limited 
to this approximation. 

In the continuum limit and with this assumption, sum- 
mations over k can be converted to integrals as follows 



/vv 
deJ2S{s-e,)\V,\' 
k 

/W p pW 

g{e)V^{e)de =- de. (13) 
'W J-w 



A. Fano-Anderson model 

Without Coulomb interaction, J7 = 0, the Anderson 
model reduces to the Fano-Anderson, or single resonant 
level, modeli^>2£. In this case the two spin channels be- 
come independent, and we need only consider the Hamil- 
tonian for a single spin species: 



H 



dU 



k 



k^k ^k 



Vk*cld), (14) 
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with omitted spin index. 

The Fano-Anderson model permits an exact sohition 
by the Green's function method^S. For a constant tim- 
neUing rate F, and provided that dot level Sd lies well 
within the conduction band^S, the ground-state expecta- 
tion value of the dot occupation number is given by 



where A{uj) = 
limit. W ^ oo 



F 



duj 



^ J-w (cj-ed- AH)'+F2' 



(15) 



^In 



we obtain 



In the infinite bandwidth 



/ X 1 1 

[nd) = — arctan 

2 TT 



if)- 



(16) 



and from the second 
5E = Ed + T 

+F 



w 



de 1 



TT SE — e 
" de 1 



w 



TT 5E + e-2£d~U' 



(22) 



We will discuss the nature of these solutions later, but 
let us note here that Eq. (j21|) from the first ansatz is [/- 
independent, and that, since these results were obtained 
from a (true) variational principle, they provide upper 
bounds for the ground-state energy of the system. 



III. APPLICATION OF CCM 



B. Self-consistent perturbation theory 

It will be instructive to compare our results for the 
Anderson model with some results already known in the 
literaturei^i^. We will compare with self-consistent per- 
turbation (SCP) results as these bear a close resemblance 
to the lowest-order CCM solutions. In Refi^, two differ- 
ent ansatz schemes are posited for ground state eigen-ket 
of the system: 



^-(11) 




IF), 



(17) 




E 



7,d|c,T ) d\ |F) (18) 



where 



IF) - E 



|0) 



(19) 



qS^kpCr 



is the unperturbed filled Fermi sea, ai.n are normaliza- 
tion factors, and rjp and are variational parameters. 
An expression for the ground-state energy of the An- 
derson model can then be obtained with the variational 
method which is here equivalent to the Brillouin-Wigner 
perturbation method^. It is convenient to present these 
results in terms of 



5E = En 



Ev 



(20) 



the difference between the actual ground-state energy. Eg 
and the energy of the unperturbed filled Fermi sea, E-p. 

For each ansatz, this procedure results in an equation 
for 5E which must be solved self-consistently. From the 
first SCP ansatz, we obtain 



5E^2Y^ 



SE 



= 2F 



' de 

lY TT 6E 



1 



ed 



-,(21) 



The application of the method starts with choosing an 
appropriate reference state 1$) and correlation operator 
S. Since for large positive ed, the impurity level has weak 
correlation with the Fermi sea, one obvious choice for a 
reference state is the uncoupled Fermi sea with empty 
impurity level 



|F)= E 



|0) 



(23) 



q^kp 



We will denote this choice as reference state I. With the 
impurity level slightly below the Fermi level, we might ex- 
pect the ground-state of the system to be approximately 
given by a filled Fermi sea plus a singly-occupied impu- 
rity. We will therefore also consider an application of the 
CCM based on the second reference state (II) 



|0) 



(24) 



q^kp,{T 



Note that these two reference state are also the refer- 
ence states for the SCP ansatz schemes in Eq. (|T7)) and 
Eq. (dH)^. In the following, we will consider each refer- 
ence state in turn, giving the correlation operators and 
deriving the CCM equations for each. 



A. Reference state I 

The correlation operators S consist of terms with oper- 
ators promoting electrons in reference states to vacant or- 
bitals. Assuming particle-number conservation, the vari- 
ous terms of S can be classified with the number of holes 
created in the Fermi sea, n, and the number of electrons 
created on the dot, m (the remaining n — m electrons are 
promoted to states in the continuum above the Fermi 
level). For reference state I with Fermi level at ep = 0, 
we may write the complete CCM correlation operator as 



N min(n,2) 



E E ^1 

n—l m— 



(n,m) 



(25) 
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with 

^(n,0) 



o(n,0) _ \p (Ln,0) TT t „ 



(26) 



^1 



(n,l) 



^1 



(n,2) 



E 



J=2 



i=3 



(28) 



The summations are performed over all relevant combi- 
nation of indices with p > kp,q < kp, with kp being the 
Fermi wave-number. 

An exact description of the interacting U ^ model 
would require that we keep all the above terms in the 
correlation operator. This is impractical, however, and 
to make progress we must truncate S in some fashion. 

The simplest nontrivial truncation scheme is to keep 
only the single term. 



Si — S'j 



(1,1) 



(29) 



We will refer to this scheme as the 5*1 or SUB-1 approxi- 
mation, and for simplicity, relabel the relevant coefficient 
as 



(1,1,1) _ Ji) 



(30) 



With this truncation of the S'-operator, the CCM expres- 
sion for the ground-state energy of the Anderson model 
reads 

E, = (0|e-^^7Je^^ |0) = 2 ^ + ^ y,;^^ . (31) 



From Eq. ([5]), the coefficient is determined by the 
equation: 

V, + [ea - e,) - ^ V^.s'^^i = 0, (32) 

q'<kF 

which is a system of A^e coupled quadratic equations. 
This should be compared with the work in Refi^ on 
the periodic Anderson model, in which Ng independent 
quadratic equations for each parameter were described. 

Under the assumption of no external magnetic field, 
the coupling constants Vk are spin symmetric such that 
the correlation coefficients can be written 



Using these new parameters Sq^' , the ground state energy 
correction from Eq. (|3ip becomes 



(33) 



(34) 



With spin symmetry, rearranging Eq. ([32]) , we obtain an 
expression of s^q^ in terms of 5E, 



V„ 



(35) 



Multiplying both sides of this equation by V* and sum- 
ming over g, we obtain a self-consistent equation for 5E: 



5E = 2 



g{e)V^{e) 



For a constant F = I{g{£)V'^ {e) , 
r° de 



6E ^ 2T 



w 



IT 5E/2 + £-£d 



(36) 



(37) 



In this form, this is like the SCP result, Eq. pTjl . except 
that the SCP result does not have a factor one-half in the 
denominator on the right hand side. This factor results 
from double dot occupancy, as will be discussed in the 
results section. 

Using Eq. ([9]), the bra-state coefficient s"q^ is obtained 
within the SUB-1 approximation from 

q' ^ 

This can be solved for s^}'' by first rearranging such that 
we have 



(1)2;(1) 



;(1) _ « V ^1' 9' 9 



SE 
2 



(39) 



and then substituting this equation into itself and iterat- 
ing. Denoting 



we find 



g(i) ^ 



V* 



- [Sd + £q) ' 



(40) 



(41) 



Therefore, once Eq. (j37[) is solved for SE, we can imme- 
diately obtain the correlation coefficients from Eq. ([35|l 
and Eq. (j4T|) . The impurity occupation can then be easily 
calculated from Eq. (O as 
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(42) 



We will also consider a more advanced approximation 
based on the Scheme I reference state. Here we choose to 
keep both terms in S with a single hole in the electron sea 
and, in order to account for the effects of the Coulomb in- 
teraction on the impurity, we will also include the lowest- 
lying double-occupation term. We therefore consider the 
correlation operator 



qa- 



q _ c(ia) , 0(2,2) £.(1,0) 



(43) 
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where 5*1^'^^ promotes electrons from the Fermi sea to the 

dot, sl^'^^ induces electron-hole correlation within the 

reservoir, and 5*1^'^^ promotes two electrons of opposite 
spin to the impurity. For notational convenience, as in 
Eq. pO|) . we relabel the CCM coefficients as 



s(i) 



.(1,1,1) 
'go- 



^ „(I-2,2) 
qq' era' qq' aa' ' 



sib) 

pgcr 



s'^'K (44) 



where a ^ a' , because of the Pauli principle. Spin sym- 
metry means that, together with Eq. ((33|) . the coefficients 
can be rewritten as 



(a) 
qq' 

(&) 
pq 

(a) 



(a) 

S^'^ - 
><3T ~ 

(a) 



(a) 

pql' 



(45) 
(46) 
(47) 



With this correlation operator, we obtain the same ex- 
pression for the ground-state energy correction as before. 



(48) 



Using Eq. ([5]), we find that the CCM expression for 
the ground-state occupation number operator, = 



- \^ zW + S''"' s^"' 



qq 



(49) 



Evaluating Eq. ([8]) for the above reference state and 
truncation scheme, we find that the ket-state CCM coef- 
ficients are determined by the following equations 



. fit''' 



ed-Sq-Y.K'^['^]+T.yp' 

q' / P 



r* „(1) I 



2ed~eq~eg,+U~2j2Vq*"s[V 

q" 

Analogously, from Eq. the bra-state parameters are 



obtained from 

s(i) Ji) 

q' 

'4n*E4'J"4'i"-4E^;'4' 



^* _ ~(i) (1 
^q *q 9' q' 



'■d — £q 



2 ' « 



* s('»)„('') 
q^q'q" 



q'q", q'q" 
5(a) ,(1) , ~(b) I T/* T/* (b) 



+4^E4?4''+E4^K;-EW 

q' P \ q' 

Vp^P + (2£d - - + [/ - 5iJ) 

^E*4"^+E^?V4^ 



(53) 

0, 

(54) 



= 0, 

Vt)-F(e,-e,)sW-F;E4'^4?- = 0-(55) 

9' 

To further simplify our equations, we scale all energies 
with the rate F, and use a bar to identify scaled quanti- 
ties, e.g. Vk = Vfe/r. We then scale the CCM coefficients 
as follows: 



VqieMe), 


4;^) = y*qieM^^ 




(a) 
%(e')q(e) 


= yq{6')yqie)^ai£',£), 


(56) 


q{e')q(e) 


= Vl{e')y*qie)^^(^''^)' 




p(p)1(^) 


= vl{p)Vq{6)£.b{p,e), 




~(b) 

'^p{p)q{e) 


= yp{p)y*q{e)ib{P,£), 


(57) 



to obtain a much simplified set of expressions. Proceed- 
ing to the integral representation (Eq. (|13p ) for the sums, 
the ground-state energy correction SE = 6E/T reads 



SE='- 

TT 



(58) 



w 



Similarly, the dot-occupation mmiber (Eq. (|49)) ) becomes 

{nd) = 2 /" -US)^i{e) (59) 
J-W I" 

+4 / — / — ^a(e',£)ea(£',e). 

Finally, the integral form of the CCM Scheme I equation 
system Eq. (f50|) - Eq. ([55|) becomes 



Ci(£) 



£d - e - 



d/9'6(p',£) 
+- r de"Ue'\e) + l 



[2ed -e~e' + U~SE] ^e, e') + C/a(£)a(£') 
^ rf£"Ki(£)^a(£',e")+6(£')^a(e",£)] 



a(£) 



1--/ d£"6(p,e") 



+ (p-£)6(p,£) 



0,(60) 

0,(6f) 
0,(62) 
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and 



1 - ?Jd + - — I 4i (£) 

- /■° de' 
+W _a(£-')ea(e',e) 
J~W 

de' f° de" 
-4 /_— _—Ue",e')Ue',e) 
-w J -w 

I" V J-W 7^ 

a (£2) + (2£d " £i - £2 + C/ - <^-B) (^l, £2) 



-w 



(e) + (p - e) (p, e) - / — (e') (p, s') 



de' 



0,(63) 

0,(64) 
0,(65) 



w 



for all £, £1 , £2 ^ and p > 0. 



B. Reference state II 

For the second reference state, Eq. ([M]) . the correlation 
operator can be written as 

N 

^(") = ^ + 4"'^) + sir''"'' + sir'''' 

n=l 



(66) 



with 



4-^^ - E 4"r.sn4.^.%^,-^ (67) 

c(".2) _ V- (ii,n,2),t TT t c r68) 

n 



{PQ<^}r. 



J = l 



4"'^ - E 4"r.£4^ini..-..^.' (70) 



where, as for scheme I, n labels the number of holes cre- 
ated in the Fermi sea. In this expression, the first term 
simply creates electron-hole-pairs in the leads, the second 
also adds a second electron to the impurity, the third 
term performs a spin-flip of the impurity electron, and 
the fourth creates the doubly-occupied impurity state. 

In this section we will consider the simple truncation 
scheme for the second reference state in which we trun- 
cate 



s=sr'^+sr'\ 



(71) 



such that we keep excitations of single spin-up electrons 
from the below Fermi surface to the dot (S'j^'^'') and of 



the spin-down impurity electron from the dot to the sea 
above Fermi surface (5*1^'°^). For this truncation scheme 
we relabel the parameters 



(11,1,1) ^ Ma) 



(11,1,0) _ (la) 



(72) 



With these choices, under spin symmetry, Eq. ([5]) yields 



Ma) _ 



Mb) _ 



V* 



E,'^;5(if)-f £,~£rf-[/■ 



(73) 
(74) 



Defining SEa = Y.p' Vp'^'"'^ and SEt = E,' V*,s(^''\ re- 
arrangement of Eq. (|75|) and Eq. (|7^ gives the two self- 
consistent equations for SEa and SEi,: 



SEa = r 
SE,, = r 



w 



de 



1 



TT SEa - £ + £rf 

° de 1 



w 



TT SEf, + e ~ ed - U 



The ground-state energy correction is in terms of these 
quantities: 



SE^ed + SEa + SEb 



(75) 



Like the solution with SCP in Eq. (p2|) . this SE involves 
two integral expressions. However, here the two integrals 
are contained within two separate self-consitant integral 
equations, whereas in the SCP expression Eq. ([22), there 
is just a single equation for SE. 



IV. NON-INTERACTING MODEL 

The non- interacting model (Eq. 1131)) is bi- linear in 
fermionic operators and thus can be exactly described 
by a CCM wavefunction with a bi-linear correlation op- 
erator. The exact correlation from reference state I for 
the Fano- Anderson model is therefore 



sr'Us\ 



(1,0) 



(76) 



In this case, the exact CCM equations are as for the 
full Anderson model in scheme I, except that C/ = 0, 
^99' ~ Oi and the spin summation is to be suppressed. 
From Eq. ([55)) and Eq. ([55)) . then, the ground-state en- 
ergy difference and impurity occupancy are 



r° de 
SE = I ^Ci(£-); 



w 



{rid 



° de~ 



(77) 
(78) 
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Under transform Eq. p3|) and scaling Eq. (|56p . the two 
systems Eq. (pO)) - Eq. (pS]) take the form 



W 7- 



1-/ — 



0,(79) 
0,(80) 



l + 6(^)(£d-e-'5S)- (nd) 



w 



dp- 



1 



" — 6(p-e-') 



6(e) 



-6(e')6(p,e') + (p-e)6(p,e) 



0,(81) 
0,(82) 



— 6(p,e) 

TT 

-W ^ 

for all £ < £f = and p > = 0. 

The number of integral equations of this system can be 
reduced by expressing ^f, in terms of ^ , and in terms of 
6 and ^i. This can be done by first rearranging Eq. 180]) 



into 



6(p, 



ei(g) 



6(p,e') 



(83) 



and substituting this back in Eq. (|80|) again to obtain 



^bip,£) 



X 1- 



■fy TT £' - p 



— 6(p,e-") 



w 



TT 



;(84) 



by repeating these substitutions and rearrangements for 
infinitely many times, we arrive at 



6(p,e) 



£ - P 1 + X(p) ' 



(85) 



where 



X(p):=x[a](p):= (86) 

^y TT £ — p 



Similarly, 
6(P,£) = 



6(e) 



1 d£-'6(£-') 



1 + X(p) TT £' - p 



(87) 



Finally, with Eq. ^ and Eq. dST]), system Eq. ^ 
Eq. dM]) is reduced to 



a(e-) 



Ed ~ s ~ 6E 



dp 



1 



/o ^ (£-p)(i + x[a](p)) 

- /■" d£' 

1 + 6 (e) (£-d ~e-5E)- —a (£-')6 (e') 

^^dp 1 1 



a(e') 



(88) 
1 = 0, 



7^ 1 + X(p) 
1 



(89) 



X Ui(^") 



de' Ue') 

1 + X(P) J-W TT £' - p 




N=25 

^N=200 
^SUB-1 
* * exact 



-2 2 

£d/r 

FIG. 1: Main panel: Ground-state impurity occupation 
number of the Fano- Anderson model as a function of the dot 
level position Ed with Fermi level at ef = 0. Plotted are 
the exact Green's function result (Eq. (|15|l 'l and results for 
the CCM Si- and S2-correlation operators. The S2 results of 
Eq. (|78|l are labelled with A'^, the number of nodes used in 
discretizing equation system Eq. ((SH) - Eq. ((SUl. The SUB-1 
curve (5*1 result) is the analytic solution of Eq. (|90|l . Inset: 
Investigation of numerical convergence of CCM-5'2 scheme. 
Here we plot (5n^^' , which is the absolute error between (rid)^ 



and the exact (ud) , 5njj = '^(nd) — (rid) |, as a 

function of N. Results are shown for two impurity-level posi- 
tions: Sd/T = and -2, with N in the range 20 — 100. These 
results show the good convergence of the CCM S2 scheme to 
the exact result. 



Equation ([89)) is a linear integral equation for S^i, so 
once a is known, ^ can be obtained directly. The equa- 
tion for a (Eq- (|88p ). however, is non-linear and must be 
solved using a numerical scheme for solving non-linear in- 
tegral equations. We use the A^th order Legendre-Gauss 
Quadrature rule^ to discretize integrals in Eq. (|88p . thus 
reducing the problem to that of solving N coupled non- 
linear algebraic equations. As exact results are available 
for the Fano- Anderson model, we are able to compare 
with these results and determine the accuracy of our nu- 
merical procedure. 

Figure [T] plots the ground-state impurity occupancy 
(ud) of the Fano- Anderson model against dot-level posi- 
tion Sd using both Si and 5*2 CCM correlation operators. 
Also plotted is the exact result of Eq. (|15p . 

The occupancy in the S'l-approximation is given by 
Eq. ([1^ but with no factor 2 in front: 



{rid) = ■ 



:(1)S(1) 



(90) 



with SE and the correlation coefficients given by their 
previous expressions (e.g. Eq. (|37p for SE) but with de- 
nominators 5E/2 + e — Ed replaced by SE + e ~ Ed to ac- 
count for the fact that we only have one spin species here. 
The S'2-solutions are the numerical solutions of the inte- 
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gral equations Eq. ((88|) - Eq. ((89|) , with different numbers 
of discretization nodes, N, in the range —4 < e^/r < 4. 

These numerical solutions are denoted {ud}''^''- In the in- 
set of this figure, we plot details of the convergence of the 
CCM 5*2 calculation to the exact result as a function of N 
{N = 20, 30, • • • , 100). We plot results for ea/T = 0, -2; 
CCM results for values of Ed away from zero always show 
less error than at Sd = 0. In this inset, wc plot the abso- 



„ g.-l □ {3 

-^^o-0- e -o- e -e- e so- 



lute difference S'l 



AN) 



between the numerical results and 



the exact solution, 



dn 



{N} 



(rid) 



(N) 



(exact) 



(91) 



These results illustrate convergence of the numerical 
results of Eq. ([75]) to the exact solution, as the num- 
ber of discrete nodes increases. From the inset, we can 
determine that for the worst case at Sd around 0, the 
absolute error is about 0.5% for N = 100, and the rela- 
tion between the absolute error and N is approximately 
error oc 1/7V^. This afhrms the equivalence between 
Eq. (1751) and Eq. (HH), and hence the fact that CCM 
is able to well reproduce the exact dot-occupation num- 
ber for the Fano-Anderson model. In the main figure, 
the SUB-1 curve is significantly different from the oth- 
ers, showing the important role of the particlc-hole cor- 
relation sl^'^\ without which the occupation number be- 
comes unsymmetric around the origin, and the dot level 
becomes narrower. 



V. RESULTS FOR FULL ANDERSON MODEL 

The presence of a finite- J7 interaction to the full Ander- 
son Hamiltonian. makes it impossible to be solved with 
a small number of correlation terms in 5*, as was the case 
for the Fano-Anderson model. In this section, we present 
our approximate CCM results for the ground-state en- 
ergy difference of the full Anderson model, and make the 
comparison with results of the SCP method. 

For our CCM calculation from reference state I we will 

(2 2)' 

keep the double occupation term Sf ' in the correlation 
operator but drop S'}^'"'' for simplicity. The correlation 



operator therefore reads 



S2a — 



(1.1) 



5i 



(2,2) 



(92) 



Figure [2] shows the ground-state energy correction 
5E, defined in Eq. ^U\\ as calculated from the various 
methods under consideration here. The results CCM I 
Si, SCP L CCM 11, and SCP II are all analytic solu- 
tions: CCM 5*1 is the self-consistent solution of equation 
Eq. (1271); SCP I and SCP II arc solutions of Eq. (gH) and 
Eq. dm respectively; and CCM II is solution of Eq. ([75l) . 
The result CCM I S2a, on the other hand, is a numeri- 
cal solution of equation Eq. (|58p with a finite number of 
discretization points = 41. From the Fano-Anderson 
calculation, we expect this solution to have an estimated 
error of 1% of F. The reason why fewer discretization 
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FIG. 2: The groundstate energy difference of the Anderson 
model with U/T = 6, W/F = 60 Ref,^. The solid line is 
numerical result of the CCM with reference State I and cor- 
relation operator S2a- N — 41 discretization nodes were used, 
with an estimated numerical error within < 1%. The dashed 
line with triangles and circles are the CCM Si solution of 
Eq. (133 (reference state I) and and CCM II Eq. ^ (ref- 
erence states II). The remaining curves are results from self- 
consistent perturbation method, with both ansatz I Eq. (|2H) 
(diamonds) and II Eq. (|22p (squares). The occupation of the 
dot at any given values of ed is roughly given by the slope of 
SE curve at that point. One therefore expects that the exact 
result has slope ~ in region > 0, 1 in — i7 < Ed < 0, 
and 2 for Sd < -U. Solutions CCM ISaa, CCM II and SCP II 
are the only approximate solutions which capture this physics 
adequately. In addition, the SCP I and SCP II results provide 
exact upper bounds for SE and we see that CCM I S2a curve 
follows this bound closely thoughout the range and thus is 
seen to be the best result in both aspects, in all regions. 



steps are used here as compared with the Fano-Anderson 
model is that for the Anderson model, the coefficient s^"^ 
has two energy indices, whereas the final equation for the 
Fano-Anderson model only involves single-index. 

The results of FIG. HI can be understood physically as 
follows. For a vanishingly small coupling between reser- 
voir and impurity (F — > 0), we can expect SE to consist 
of three straight lines: for e > 0, we expect there to be 
no electron on the dot such that 5E ~ with slope zero; 
for — [/ < Ed < 0, one electron from the Fermi sea will fill 
the impurity, yielding SE ^ ed (slope 1); and finally, for 
Ed < —U , two electrons with opposite spins will overcome 
Coulomb interaction and occupy the impurity, such that 
5E ~ 2ed + U (slope 2). For finite F, we expect the ex- 
act solution of SE to broadly exhibit the above features, 
with transitions between them over an energy scale of or- 
der ~ r. In the figure, CCM I S2a, CCM II, and SCP II 
curves all demonstrate three regions where the solution 
has gradients 0, 1, and 2. In contrast, the CCM I re- 
sult only has and 2 gradients, whereas SCP I has only 
and 1 gradients. This difference in the negative Ed regime 
is a manifestation of the exponential structure of a CCM 
ground state ansatz — although there is no double occu- 



9 



pation correlation in 5*1 , having Sf in its expansion e'^^ , 
it can still produce a doubly occupied dot. 

Moreover, since SCP method is equivalent to a varia- 
tional calculation, the SCP results give upper bounds for 
the exact value of SE. In the region £d > 0, SCP I gives a 
lowest upper bound, whereas for Ed < 0, it is SCP II that 
gives the lower one — neither solution provides a consis- 
tent upper bound for the whole range of e^. However, as 
is clear from Eq. ([2]) the CCM II solution follows rather 
closely the best upper bound from both SCP calculations 
across the complete range. This distinguishes it as, of all 
the approximate methods discussed here, as producing 
the best result across the whole parameter space of the 
model. Note however, that the CCM solution does not 
provide an upper bound. 

VI. CONCLUSIONS 

Our calculations demonstrate the applicability of the 
CCM to quantum impurity models. Comparisons with 
exact results for the non-interacting Fano- Anderson 
model demonstrate that the CCM can determine ground- 
state properties for such models to a high degree of ac- 



curacy. 

In typical applications, the (truncated) CCM normally 
works well when the exact ground-state of the system re- 
sembles the reference state to some extend. Here, how- 
ever, the method is able to produce good results for the 
(interacting) Anderson model in all regimes; even the 
simplest truncation scheme is able to describe double oc- 
cupancy of the impurity. 

In our calculations, the main technical difficulty is with 
solving integral equations, while higher-order CCM coef- 
fcients will have more wave-number indices (j> and q). If 
we have N discretization notes for each integral, n in- 
dices require N"- data points per coefficient. This prob- 
lem can be reduced by resummation of electron-hole con- 
tributions as is done here for the Fano- Anderson model. 

Future work includes the extension of our calcualtions 
to multiple impurity problems, such as the periodic An- 
derson model^, and to include higher-order terms in the 
correlation operator to fully cover the Kondo effect. 
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